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ABSTRACT 

The ordering dynamics of the Higgs field is studied, using techniques inspired by 
the study of phase ordering in condensed matter physics, as a first step to under- 
standing the evolution of cosmic structure through the formation of topological de- 
fects in the early universe. The common feature of these different physical processes 
is scaling. A fully analytical approximate scheme — the linear-gaussian approach 
— is proposed to evaluate 1-point, 2-point, etc. scaling functions for the ordering 
dynamics of the 0(?2) -symmetric Higgs-field models. 
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I. INTRODUCTION 

When a system is rapidly quenched from a disordered phase of high symmetry to a 
multi-phase region of lower symmetry it undergoes a spontaneous symmetry break- 
ing (SSB) phase transition. During this transition the system develops a spatial 
structure of randomly distributed domains which grow with time. This phase order- 
ing process has been extensively studied in the context of condensed matter systems 
especially those with a non-conserved order parameter (model A) 0], described 
by the time-dependent Ginzburg-Landau (TDGL) equation. There is much evidence 
that in the late stages of growth these systems enter a scaling regime 0, in which 
the two-point correlation function has the scaling form 



where is the vector order parameter field, L{t) is the characteristic length scale 
at time t after the quench, / is a scaling function, and angled brackets indicate an 
average over initial conditions (and thermal noise, if present). 

A similar kind of phase-ordering phenomenon is believed to have occurred in the 
early universe. While the big-bang theory has been widely verified by observations 
(confirming that the universe began in a very hot, dense state and has expanded 
and cooled down ever since 0]), the origin of cosmic structure remains unexplained. 
According to the isotropy of the microwave background radiation (left over from 
the early matter-radiation decoupling transition) the early matter distribution was 
very uniform. How did the universe evolve from its primordial smooth state to its 
present state of lumpiness, where matter concentrates in galaxies and galaxy clusters 
1^ forming a very-large scale structure? It is believed that tiny large-scale density 
fluctuations, present at decoupling, could, if strong enough, have resisted the overall 
expansion and grown under gravitational collapse. Matter in the overly dense regions 
of space would have clumped together to produce general structure. What was the 
origin of these small fluctuations, however, and how could they have generated the 
kind of large-scale structure we see today? Based on a process central to unified 
theories of particle physics — that as the early universe cooled down a hypothetical 
field, the Higgs field, underwent a SSB transition — it has been suggested that the 
consequent field ordering and defect formation could have provided the mechanism 
to generate structure. Field defects would form unavoidably, because 'vacuum' 
configurations above the horizon scale are uncorrelated. Since the defects carry 
energy they could provide the fiuctuations around which matter would aggregate 



The purpose of this paper, is to use some of the techniques developed in the 
framework of 'model A' dynamics (i.e. the TDGL equation) to study the Higgs 




(1) 
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model ordering kinetics. This problem is technically more difficult than model A 
because the equation of motion describes a damped wave propagation rather than 
a purely dissipative process. However, these non-conserved field ordering processes 
are likely to exhibit similarities at late-times, where a scaling regime is expected to 
occur 1^, |10|, |Tl]]. A difference, though, is that here the characteristic length scale 
grows linearly with time, L(t) ~ t, while L(t) ~ \/i for model A. 

While domain growth phenomena, governed by the kinetics of topological de- 
fects, have been fairly well understood within model A dynamics, a first principles 
calculation of the scaling functions has proved to be a most difficult task, and various 
closed-approximation schemes to evaluate the scaling function f{x,q) of equation 
(|l|) have been proposed in the past few years [|^, The key technique, exploited 
by several authors [|l2|-||l6l, is to to introduce a mapping 0(r,t) = 0(m(r, t)) be- 
tween the order parameter field and an auxiliary field m(r, t) which has, near a 
defect, the physical interpretation of a position vector relative to the defect. With 
this new variable, the problem of describing the field at each instant of time is 
transformed into a problem of describing the evolution and statistics of the defect 
network. This approach enables the use of a physically plausible and mathemat- 
ically convenient gaussian distribution for rh. Such a distribution is unacceptable 
for the order parameter itself, since this is effectively discontinuous at the domain 
size scale. The application of this sort of approach to the scalar-field model A has 
recently received a critical review by Yeung et al. [0]. Methods based on a de- 
scription of the wall dynamics lead to an approximate /m ear equation for m{r,t), 
or for its correlator (m(x, i(:i)m(x + r, t2))- A different and promising approach, due 
to Mazenko [Q, aims at deriving a closed non-/mear equation for C(r, ti, ^2), built 
on the equation of motion for the scalar-field model A, and the assumption that the 
field m is gaussian distributed at all times. It has the virtue of yielding results with 
a non-trivial dependence on the spatial dimension d and it is also easily extensible 
to 0(ra)-component systems. Despite the uncontrolled nature of the gaussian as- 
sumption these approaches have been shown to give good results, displaying most of 



the expected physical features |jTj, |T6|]. For the nonconserved dissipative dynamics 
of model A, it has been argued that the gaussian approximation becomes exact in 
the limit of large spatial dimension d, while for fixed d it provides the starting point 



for a systematic treatment [0. It is also correct for any d in the limit of large 
n. For the Higgs field model considered here, the gaussian approximation is again 
exact for large n, but the large-c? limit does not seem to be simple. Nevertheless, 
by incorporating topological defects in a natural way, the gaussian field approach 
provides the simplest non-trivial approximation scheme for the dynamics of phase 
ordering. 
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In section 3 we attempt to apply the Mazenko approach to the 0(n)-field Higgs 
model. The late-time pair correlation function is then given by the Bray-Puri-Toyoki 
(BPT) function (pOD ||15|, as a function C{'~f,n) of the normalized correlator of fh, 
7(r, ^1,^2), which obeys an approximate non-linear equation. The BPT function 
embodies the asymptotic defect structure, while 7(r, ^1,^2) describes the dynami- 
cal dependence of C(r, ^1,^2)- The mapping used by Mazenko []T3|, |T6[, however, 
restricts the field to evolve within the bound |0| < 1, which is incompatible with 
the oscillatory bulk relaxation of the Higgs field, and leads to an inconsistent ap- 
proach. The difficulties with this approach, however, motivate our next attempt to 
tackle the problem. In section 4 we consistently eliminate the field bulk oscillations 
by restricting the asymptotic field dynamics to the 'vacuum manifold'. Extending 
Mazenko's gaussian approach to the Non-Linear Sigma Model (NLSM), using the 
unit vector mapping (p = m/\m\, the pair correlation function is still given by the 
BPT function, but now 7 obeys a different approximate equation. Rather then 
solving numerically this equation for 7, which is rather complicated, in section 5 
we propose a fully analytical approximate scheme to evaluate C. This amounts to 
replacing 7 by its large-n limit in the argument of the BPT function, but keeping 
the remaining n-dependence unchanged. As n — 00 the equation for 7 becomes 
linear and exactly solvable 0, so we may call this scheme the Linear- Gaussian (LG) 
approach. Although we cannot use the NLSM for a scalar field, the approach still 
holds for this case if one takes n = 1 in the BPT function. In section 6 the LG 
approach will be generalized to evaluate other kinds of scaling functions, such as the 
average of the energy density and its correlation function. 

II. THE HIGGS FIELD MODEL 

In this section we briefiy review basic notions about the cosmological background 
model. The Higgs field model is presented and its dynamics are discussed. 

As is usual practice, we shall consider a flat expanding universe as the model for 
the bulk cosmological background @, 0. In this case the local curvature is zero and 
the metric is space- independent, given by 

ds^ = cUt^ - a{t)V = a{r]y (^c^ df]^ - dr^) , (2) 

where: t and r are comoving coordinates (i.e. the reference frame is moving with 
the cosmic fiow); a(t), or a{ri), is the space expansion factor; 77, the conformal time, 
defined by 

drj = dt/a{t) , (3) 

plays the role of 'real time' in a static universe: the horizon of an 'event' after a 
time t — its maximum range of infiuence after time t — is given by h{t) = Jq dr = 
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/q c dt/a{t) = c Jq drj = c rfit). In a flat universe, the function 

varies slowly with time from a = 2 (radiation era) to a = 4 (matter era). Away 
from the matter-radiation decoupling transition a can be regarded as a constant 
and the expansion factor is given by a power-law a ~ ^ 7^°/2 ^ 

In the early radiation dominated era the energy was dominated by relativistic 
particles (with equation of state p = p/3), yielding a ~ t^^'^ ~ r] and p ~ 
~ t^^. Here p and p are the uniform background pressure and energy density. Once 
the universe cooled down and matter decoupled from radiation this became the 
dominant source of gravitation (with negligible pressure p ^ p), yielding a ~ 

^2/3 

~ Iff and p ~ ~ t"^ in the matter dominated era. As matter became transparent 
to radiation, the matter perturbations started to grow. 

A simple class of SSB theories is provided by the (real) n-component Higgs field 
models, where a 'global' 0{n) symmetry is broken [^, |^, |Tl[]. These theories include 
several cases where topological defects form: domain walls {n = 1), global strings 
(n = 2), global monopoles {n = 3) and global textures (n = 4), which are of potential 
interest as a mechanism to generate cosmic structure. 

The dynamics of the Higgs field 0(r, t) = (0^, 0") in an expanding universe is 
derived from the Lagrangian density |^ 

C{r,t) = i {d^/dtf - \ (V0/a(t))' - r(0) , (5) 

where V is with respect to comoving coordinates, and 1^(0) is a generalized 'double- 
well' potential with an 0{n) symmetric 'vacuum manifold' where 0^ = 1. Minimizing 
the action S = J dt J d^r a{t)^ C{r,t) (where dtd^ra{t)^ is the covariant 4- volume el- 
ement) with respect to variations of 0, and using conformal time, yields the equation 
of motion 

7r^ + -7r = ^'^-«^^' 6 

or]-' r] or] 90 

— * 

a wave equation with a damped 'friction force' {a/r]){d(p/dri), which mimics expan- 
sion in the comoving frame (and destroys the Lorentz invariance), and a non-linear 
force dV/d(f) which drives the field to the 'vacuum manifold'. The initial conditions 
— corresponding to a disordered state before the SSB phase transition — shall be 
discussed in the appendix, where we present the solution of (|) in the limit n ^ 00. 
The price for using conformal time is to have an effective potential in equation (^ 
with a time-dependent amplitude. The particular form of the potential, however, 
should not affect in any essential way the late-time dynamics and scaling properties. 
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We expect, for instance, the main effect of 0(77)^ to be a decrease by a factor 1/a 
in the comoving size of the defects core, which simply speeds up the system entry 
into the scahng regime. To simphfy the subsequent discussion we shall from now on 
discard the factor in the equation of motion (^. We will not really need that for 
computational purposes, as we shall be using the NLSM. 

Taking conformal time on the same footing as real time, equation can be 
viewed as a 'general-relativistic analogue' of the TDGL equation, describing the 
dynamics of non-conserved systems. The Higgs Hamiltonian density corresponding 

to (D 

^(r, t) = ^ { ^ (d^/dri) ' + ^ ( V0) ' + V{$) } , (7) 

is (apart from 1/a^) similar to that of a static (Minkowski) universe, and compared 
to the TDGL model has an extra 'centripetal' term [dip/drff'. For a vector field in 
the 'vacuum manifold' it leads to an energy density which decays (due to expansion 
and dissipation) like the background, p ~ Therefore, the Higgs field yields 

density fluctuations of constant amplitude 5p/p = {p — p)/p^ const which, through 
Einstein's equations, provide a source for perturbations in the matter distribution. 
Assuming the existence of a late-time scaling regime (which has been confirmed 



by numerical simulations |T0|, [T^]), the dimensional analysis of (|g) leads to a 



characteristic scale growing with the horizon 

Hv) - crj, (8) 

implying that the field defects move with relativistic speed. We therefore expect 
the pair correlation function (|l|) to take the asymptotic scaling form C{r,rii,ri2) = 
f{x,q), with scaling variables x = r/rji and q = ri2/rii, where r = |ri — is the 
distance between the two points. 

Causality constrains the field correlations after the SSB transition. Two field 
configurations at times rji and 772 can only be causally correlated if their distance 
r is below the sum of their horizons crji and 0772 (i.e. if the horizons intersect). 
Therefore, the condition for C{r,rii,r]2) 7^ is (taking c = 1) 

r < V1 + V2 ■ (9) 

If one of the horizons contains the other configuration [rji or 772 > r) the correlations 
are 'direct'. Otherwise, 'indirect' correlations can occur through common causal 
correlations with intermediate points in the region of intersection of the horizons. 

Unlike purely relaxing systems, the wave nature of the Higgs dynamics forces 
the late-time saturating field not to satisfy |0| < 1 even if its initial condition does. 
To see how the field tends to its 'vacuum manifold' we linearize equation @ as 
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approaches a given 'vacuum' state 0o. Considering a single-domain region where (j) 
can be taken as uniform, and noticing that the only restoring force is parallel to 0o 
(normal to the manifold) due to the symmetry of the 'vacuum', we find, at late-times 

0(77) -00 ^ 1 7-T {ci(a) cos{Ar]) + 02(0) sm{Ar])} , (10) 

a{7]) 

where A = [d^V/di$^f)^, and Ci(a) and 02(0) are arbitrary constants. For a 
scalar field 0(77) ■ 0o is replaced by 10(^7) |. We conclude that the Higgs field saturates 
with damped oscillations. 

III. GAUSSIAN THEORY FOR A 'SOFT' FIELD 

In this section we apply to the Higgs model the gaussian approach proposed by 
Mazenko [Q for the TDGL dynamics. Although the approach, which is based on 
an unphysical mapping for the Higgs dynamics, leads to an inconsistent theory, it 
will motivate the implementation of a gaussian approach for a unit vector field in 
section 4. 

To derive an equation for the pair correlation function (|l]), we multiply the 
equation of motion (|^), evaluated at point (1) = (ri,r7i), by 0(2) = 0(r2,?72) and 
average over the ensemble of initial conditions, yielding the exact equation 

rv 

C(i, 2) + - C{1, 2) = C(l, 2) + ^(1, 2) , (11) 
where the driving force, or non-linear (NL) term, is 



and C(l,2) = dC{l,2)/dr]i = ^0(1) -0(2)^, etc. To transform (|TTD into a closed 
equation we need to write the NL term as some approximate non-linear function 
of C(l,2). A key idea, exploited by several authors within the TDGL dynamics 
[p!2| , |I4|, [T^, is to employ a non-linear mapping between the order parameter 0(r, r/) 
and an auxiliary 'smooth' field m(r,?7). This new variable describes the late-time 
defect network structure, which will have formed at the late stages of field ordering, 
and allows for the approximation to be implemented. The most obvious way to define 
the function 0(m), is to follow Mazenko's suggestion of using the equilibrium 
profile equation of an isolated defect (in a comoving frame) 

= 9^/90, (13) 

— * — # 

with boundary conditions 0(0) = and 0(m) m/\m\ as (|m| — * 00), and where 
Vm is the gradient with respect to rh. Close enough to a defect (i.e. for |m| <^ L{ri), 
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where the field is unaffected by neighboring defects) m(r, rj) can be identified as 
the CO moving position vector of point r from the (nearest part of the) defect. This 
picture requires, of course, that n < d. With ([T3| ) the magnitude of (p{m) is a 
monotonically increasing function of the magnitude |m|, approaching for large |m| 
the 'attractor' value 1 imposed by the potential. For a scalar field, the function 
(j){m) has a typical sigmoid form. 

The mapping (|1^) restricts the field magnitude to be |0| < 1. This is appropriate 
for diffusion fields evolving from a disordered state, but is physically incorrect for the 
Higgs field dynamics, where the system self-organizes oscillating about the 'vacuum' 
states, as shown by ([ToD . While we can prove that the use of (p!3|) leads to an 
inconsistent theory [T^ , it seems unlikely that an adequate one-to-one mapping could 
be defined for this problem. In section 4 we shall overcome this technical difficulty 
by restricting the field dynamics to the 'vacuum manifold', (f)^ = 1. Meanwhile, 
for completeness we will pursue this approach a little further using ([l3|) to derive a 
closed equation for C(l, 2). 

Following Mazenko fl^, we now assume that m{r,7]) is a gaussian random field 
(with zero mean) at all times, described by the pair distribution function 

,™ , ^ /m(l)2 m(2)2 27m(l)-m(2) , , , , 
P(m(l), m(2)) = AT" exp <^ -— + - \ ' ^ ^ } (14) 
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7^)^50(1) So{2) ^Soil)Soi2) 



So{l) = {m{lf) , Co(l,2) = (m(l)m(2)) , (15) 
7(1,2) = , (16) 



where = (^n^J (1 — 'j'^)So{l)So{2)^ , and m(l) and m(2) are the same (arbi- 
trary) component of m(l) and m(2). All the averages over the ensemble of initial 
configurations are replaced by gaussian averages on m, and can be evaluated as 
functions of the second moments 6*0(1), 5*0(2), and 7(1,2). However, from (|^) and 
the mapping (|T^) , according to which fh can be identified as a position vector, we 
anticipate the asymptotic scaling form 

Soiv) = ^ - Liv)" ' (17) 

where /? is a parameter to be determined. Thus, within a gaussian approach, the 
only variable in the problem is the function 7(1, 2), which accounts for the particular 
dynamics of 0. 

The driving force (|12D in equation (|TT|) is then given, as a non-linear function of 



C(l, 2), by Mazenko's result |T|, [16 
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where = (dC/d'y) and we have used (|T7D. Note that, by use of the mapping (0), 
there is no longer any exphcit dependence on the potential V{(j)) in (|l^), though the 
relation between (p and fh depends on V. At late-times the field will be saturated 
almost everywhere except at the defect cores (whose size is much smaller than the 
domain scale), and we may, for simplicity, evaluate the gaussian averages replacing 
the profile mapping (|13|) by its discontinuous asymptotic form 0(m) = m/|m|, for 
a vector field, or 0(m) = sign(m), for a scalar field. At late times, therefore, the 
detailed form of the potential is not important (although it must, of course, have 
the 'Mexican hat' form in order to support non-trivial solutions of (|1^)). This is in 
accord with the expected 'universal' nature of the late-stage scaling behavior. 

Evaluating the pair correlation function C(l,2), using (|l^) and the mapping 
above, yields the explicit relation C = C{'~f,n), which we will call the 'BPT function' 

H. 

CW = (^-M) (19) 



|m(l)| |m(2)| 



2 



,(,2)!imillf L!i±l;,(,2)^) . (20) 



27r V2 2' 2 

where B{x,y) is the beta function and F{a, b; c; z) is the hypergeometric function. 
The substitution of (|I^) and C(7, n) into equation (|TTp yields the approximate closed 
equation for 7(1,2), which for a vector field must be regarded as the independent 
variable. In the limit n —>■ 00 the BPT function reduces to C{'~f, 00) = 7, yielding 
JF(1,2) = C(l, 2)/S'o(l), and (|TT]) becomes a linear equation. 

We now focus on the pair correlation function at equal-times [rji = = f]), 
which is of interest by itself and also yields the initial condition to solve the general 
equation |]16|. Equation ([TT|) then reads 

1 ■■ 1 ry ■ 

- C{1, 2) - C(l, 2) + -- C(l, 2) = V'Cil, 2) + ^(1, 2) , (21) 

where (7(1,2) = dC/dr], etc. The unknown quantity C(l,2) may be eliminated to 
get a third order equation in C. Then, replacing JF(1,2) by its approximate form 
([18|) , using (p!?!) , and looking for an isotropic scaling solution C(r, 77) = /(x), which 
implies 7(r, 17) = 7(0;), with x = r/r/i, leads to the equation for 7(x): 

^^i^ {7'" + SlW'D,, + y'D,,,} + l^x' + 2(rf + 1 - a)) {7" + i'D,,} 

- g(« - 2)(2a - 3) + ^{a - l){d - 1)) 7' = -(3 {(1 - «) 7 + x(l + ^D,,) 7'}(22) 

where 7' = d'j/dx, etc, -0^7(7) = C^^/C^, ^^^^(7) = C^^^/C-y, and = dC/d'y^ 
etc. The NL functions D^^{^) and D^^^^'j) are obtained from (|20|) and embody all 
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the 77,-dependence of (p2|). The boundary conditions for equation ( p2D are 7(0) = 1, 
from definition ([T5|) , 7'(0) = 0, from 7(x) = 1 — O(x^) as a; — * 0, and 7(2) = 0, from 
C(7) ~ 7 as 7 — > and the causal condition /(x) = for x > 2. We notice that 
the boundary points are both singular, which makes the numerical solution of ( p2D 
difficult. 

For a scalar field the BPT function (|20| ) can be inverted to give 7 = sin (7rC/2), 
yielding a NL term JF(1, 2) = (2/7rS'o(l)) tan (7rC(l, 2)/2). Hence we can express 
([2^ ) as an explicit non-linear equation for the scaling function C(r, rj) = f{x): 

m'" + (x^^^ + 2(. + 1 - «)) /(X)" 

- (^(« - 2)(2« - 3) + ^(a - - 1)) fix)' 
= - /5 1^ (1 - a) tan (|/(x)) + x sec^ (^/(x)) /'} . (23) 



To perform a small-x expansion of (|23D , we recall that with the mapping = 
sign(m), used to evaluate C(7), the condition /(O) = (0^) = 1 has been built into 
the theory (although in an inconsistent manner). We find that /(x) admits a series 
in odd powers of x (implying that all derivatives at x = are determined without 



recursion), giving the linear behaviour, or Porod's regime [£l 



which is a physical consequence of having 'sharp' walls at late times. To find the 
small-(2 — x) asymptotic form of /(x), we notice that as /(x) — > equation ( p3D 
becomes linear and has three independent solutions. Since the singularity at x = 2 is 
regular we try a Frobenius power series solution Ao{2 — xY {1 + J2'kLi flfc(2 — x)''), 
and find that the equation admits a leading power-decay /(x) ~ (2 — x)^ as x — > 2, 
where p can assume any of the values: p = 0, p = 1 or p = a + {d — l)/2. p = must 
be excluded as being incompatible with the boundary conditions (it would imply 
Aq = 0), and thus the solution has the general asymptotic form as x ^ 2: 

fix) ~ (2 - x)°+('^-^)/2 1^ ^ - x)} + (2 - x) {1 + 0(2 - x)} . (25) 

Since the BPT function ( pOD has the same behaviour C(7, n) ~ 7 (and D^^, D^^^ 
0) as 7 — > or n 00, to linear order in the regime x 2 and / ~ 7 equation 
(|2^ ) is n-independent and identical to its large-n limit. In the appendix we discuss 
the large-n limit of the NLSM and find that /oo(x) ~ (2 - x)"+^ as x 2, for 
d = 3. Therefore, for any value of n, and at least for short-ranged initial conditions, 
we expect the leading power-law decay 

fix) ~ (2-x)"+^ , (x^2) . (26) 
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Although we are not looking to solve equation (^), we describe how one in prin- 
ciple could do it. From (p^ and @, the boundary conditions are: /(O) = 1, 
/'(O) = —{l/Ti)\j2l3a/{a — \){d — 1) and /(2) = 0. The parameter j3 is numer- 
ically determined by imposing the coefficient of the dominant solution in (|25|) to 
vanish, Aq'\i3) = 0. In the large-n limit, where equation ( ]Tl| ) becomes linear and 
the gaussian approach is exact, (3 can be found analytically. Comparing (|1T]) with 
the linear equation (p6D, which amounts to compare the limit of (0), JFoo(l,2) = 
/?^7(l,2)/2r/2, with (r(l)) given by (0), yields 

/3oo = 2To = 3(2a+l)/2 . (27) 



In conclusion, although the mapping (|TB|) discards the field oscillations (|TD|) and leads 



to an inconsistent theory [11^1 , equation despite its intrinsic incorrectness 



bears no obvious signs of inconsistency. A Porod's regime (p^ is obtained as a 
consequence of the 'sharp' wall constraint 0(m) = sign(m) used to evaluate C(7, n). 
We have shown that the manner in which /(x) vanishes at x = 2, given by (^), is 
independent of n and exact. 

IV. GAUSSIAN THEORY FOR THE NON-LINEAR SIGMA MODEL 

In this section we study the dynamics of a vector Higgs field within the NLSM. By 
constraining the field to lie on the vacuum manifold, this model automatically avoids 
the technical difficulties associated with the asymptotic bulk oscillations noted in 
section 3. We develop a gaussian approach, analogous to that of section 3, and 
derive an approximate equation for C(l,2). 

Long after the SSB phase transition the driving potential V closely confines the 
Higgs field to the 'vacuum manifold' almost everywhere (except at the field defect 
cores). We have shown, however, that the wave nature of the dynamics leads to a 
field bulk saturation accompanied by slow decaying oscillations about the 'vacuum 
state', preventing us to define an adequate one-to-one mapping between (p and an 
auxiliary field rh. The mapping ([T3|), for instances, forces the field to obey \(f)\ < 1 at 
all times and yields an inconsistent approach. To overcome this technical problem, 
we notice that the oscillations (p!OD are unlikely to have a major effect on the late- 
time dynamics of the field defect network (and thus on the scaling properties), and 
may thus be consistently discarded by restricting the 0(n) field dynamics to the 
vacuum manifold. Replacing the vanishing driving force dV/dcf) in (^ by a non- 
linear coupling term which constrains the length of the field, the field evolution is 
now described by the non-linear sigma model (NLSM) equation 0: 

+ - ir- = VV + T(r,r/)0, (28) 



drf 1] dr] 
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where T(r, t]) is the free Lagrangian density in (|^) 

T,r.„).(g-Vv)..^^(V^T-(g)^ 

As another advantage of using the NLSM, the ordering dynamics becomes indepen- 
dent of the details of the potential V{(j)) and, in particular, the factor 0(77)^ in @ 
is suppressed. 

The exact equation for the pair correlation function C(l, 2) is still given by ( |TT|) 
C(i, 2) + - C{i, 2) = 2) + ^(1, 2) , (30) 

where, from (pSj) and (^), the NL term is now given by 

^(1,2) = (t(1)0(1)-0(2)) , (31) 

which must be replaced by some approximate non-linear function of C(l, 2) in order 
to transform ([TT|) into a closed equation. Following the strategy of section 3, we 
introduce a non-linear mapping between the order parameter 0(r, 77), which is now 
not well defined near the defects, and an auxiliary 'smooth' field m(r, ■)]). We can no 
longer define = 0(m) using the equilibrium profile equation of an isolated defect 



14| , which yields a trivial relation everywhere except at the defect cores where it is 
singular. The natural way to define the relation between the unit vector 0(r, rj) and 
7n{r,ri) amounts to replacing (ITB|) by its discontinuous asymptotic form 



Hrh) = —r . (32) 



This mapping only determines m(r, 77) up to a factor (which e.g., may be a function 
of time), and there is now no obvious physical interpretation for the new variable. 
Up to a factor, however, we may still regard m(r, rj) as a position vector (close 
enough to a defect) like in section 3. 

For mathematical convenience we assume that m(r, 77) is a gaussian random field 
(with zero mean) at all times, described by the pair distribution function (|T4|)-(|r5|). 
All the averages over the ensemble of initial configurations are replaced by gaussian 
averages on rh, and can be evaluated as functions of n, 5*0(1), 5'o(2), and 7(r, rj), the 
normalized ?7i-correlator, which contains all the dynamic dependence. 

In the same spirit which lead to expression (|l^) in section 3, using the mapping 
(|3^) and the gaussian property of rh we can shown p3| that the NL term (|3l|) is 



then given, as an approximate non-linear function of C(l, 2), by 



^(1,2) 



{(.(l)^)-((V.(l))^)}2^ + 5o(l)^i^ 
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+ S„(l)a(i.2)l{-^^;^j||^H-Q„(l,2)} (33) 
Q„(l,2) = (n-3)(n-l)((|m(l)|'|m(2)l)"') (34) 



2 



5o(l)v/^o(l)^o(2) 



2iT V2'2'2 



F -,;^;-;7(l,2)^ (35) 



Using (|T5D, ([T6|) and (pO|) the NL term (|33|) can be fully expressed in terms of 
7 and Sq. For example, the correlators (m(l)^) = Co(l,2)2_*i and ((Vm(l))^) = 
-V2Co(l, 2)2^1, the derivative dC/dSo{l) = -C^ l/2So{l), where C^{-f) = dC^/d-f, 
and similarly for the other derivatives of C(l, 2). Substituting the NL term and the 
BPT function (^) into ( pOD we get the equation for 7(1,2), which is the indepen- 
dent variable. Specializing to equal-times {rji = ri2 = t]), and looking for an isotropic 
scaling solution 7(r, 77) = 7(x), we then obtain an approximate closed equation for 
7(x), the NLSM version of (P^, the boundary conditions for which have been given 
in section 3. Even if we take Sq to be time independent, this equation will still be 



much more complicated than (^4 ). 

If fh in the NLSM is set to have the same interpretation as in section 3, and thus 
to obey ([T7|) , we may compare the NL terms (0) and ([T8|) . The NLSM gaussian 
approach generates the 'soft' field result, as long as ((Vm)^) = 1, plus additional 
terms following from the consistent use of the mapping (132|) . This differences indicate 
that the gaussian approach is not quantitatively accurate, since (H) and (pSD should 
yield equivalent asymptotic dynamics. 

V. LINEAR-GAUSSIAN APPROXIMATION 

Rather then solving the extremely complicated approximate non-linear equation for 
7(1, 2), we propose a fully analytical scheme — the 'Linear-Gaussian' (LG) approach 
— to evaluate C(l, 2), which combines a gaussian mapping for a unit vector (p with 
the large-n exact solution. 

We notice that the relation C = C{'y,n), defined by (|TU|) and given by the 
BPT function ( ppj ) for a gaussian m, accounts effectively for the presence of the 
field defects (through the orientation of = m/|m|), and also for their topological 
nature (through the ra- dependence), and so it already describes fairly well the late- 
time defect structure. Hence, the particular form of the function 7(1,2), which 
contains the dynamical dependence of C(l,2), should not be so relevant and may 
be approximated rather crudely. For simplicity, we replace 7 by 700, the exact 
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solution in the large-n limit p5[. The scaling function fLci^s^ci) = C{1,2)ig with 
n = 1, ..,4, oo and a = 2 and 4, obtained using this procedure, is plotted in figures 
1 and 2 with fixed values of q = ri2/i]i and abscissa Xs = 2r / [rji + ri2) , and in figures 
3 and 4 with fixed values of Xg and abscissa q. More details are given in section 7 
and in the figure captions. 

As n cxD, \rn\ = {Y^i'mf)^^'^ —>■ y/nSo, and it is easy to find the limit of the 
functions Qn and C(7), either from their definitions ( p4D and (|19|) or from their 
gaussian averages (pS]) and (pO[). The BPT function reduces to C(7, ra) — > 700 = Coo 
and Q„(l,2) ^ l/^o(l) v/'S'o(l)^o(2). Equations (|g)-(|31|) yield the self-consistent 
linear equation 

7oo(l, 2) + - 7oo(l, 2) = 7oo(l, 2) + (T(l)) 7oo(l, 2) , (36) 

(T(l)) = n/vl , (37) 

where the scaling form (|37D follows from a dimensional analysis of (|36|) or (^9]) (and 
from translational invariance), and the constant To is to be found self-consistently 
(see (|6^) in the appendix). The linear term J^ooi^, 2) = {T{rii)) 7oo(l, 2) is the limit 
of the previous NL term: the gaussian expression (^), or the definition (|3TD and 



(p9|), where m/ v^TSq. 

Instead of determining 7oo(l, 2) by solving the linear equation (P^ ) at equal times, 
which (like (PT|)) is third order, it is easier to calculate the correlation function of the 
exact large-n solution of the NLSM equation (PB]), which is second order. Equation 
(P^) for 7oo = Coo can be derived from the large-ra limit of ( ^5]) (just like (^Up was 
derived from (^)), so the two procedures to obtain 700 are equivalent. Turok and 
Spergel have solved the large-n NLSM in momentum space and determined the 
structure factor corresponding to a random initial field. In the appendix we present 
and Fourier-transform their result to 3-dimensional real space, yielding the scaling 
function /oo(x, q) = Coo(r, r^i, r^a) = 7oo(r, Vi, m)- 



9{l + q-x) 1 .2nh^/2 / 



a + l 



fooix,l) = ^^ I. dss{l-s')^{q'-{x-sY)^ , (38) 



where 



x = r/r]i , q = r]2/r]i , (39) 



{B,A) = {x + q,x — q) , x<l — q 

= (l,x — g) , \l — q\ < X < 1 + q 

= (1,-1) , x<q-l, (40) 

and N = 4/5, 32/63 for a = 2, 4. At equal-times, 

fooix, 1) = i /' ds s{l - s')^{l - (X - sf)^ . (41) 

JM X Jx-i 
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The small- a; behaviour of foo{x, 1) can either be obtained from the large-n limit 
of the gaussian equation for 7, e.g. (^2]) , or by expanding (^l]) as x 0. Both 
procedures yield the leading behaviour as x — > 

f^{x,l)= 1 - (5/8) In (1/x) a;2 + ■ ■ ■ a = 2 

1 - (27/16) x2 + • • • a = 4 . (42) 



Expanding the BPT function (^) as 7 — > 1 and using ( fl2[ ) yields the small- a; 
expansion for the equal-times pair correlation function within the LG approach. 
For a scalar field we have 



fLGixA)= 1- W51n(l/x)x + --- a = 2, n = l 



= l-lp7/2x + --- a = A, n=l, (43) 

and for a vector field 

fLGix,l)= l-Ai{x)x'^-\ a = 2, n>l 

= l-AaX^H a = 4, ra>l, (44) 

where Ai{x) = (5/8)(ln(l/x))2 and A2 = (27/16) ln(l/a;) for n = 2, and Ai{x) = 
(5/4) ln(l/x) and = 27/8 for n = 3. Performing a small- (1 + q — x) expansion of 
we find the leading power-law decay 



In the limit of very-different times [rji <^ 772)? we obtain the leading time-decay 



1 Q ( ^±1 3 /2 ) 

f{x,q) ^ Ux,q) ^ B {1,3/2) + ^^^/^'^') ' - 00) . (46) 

By the same arguments discussed in section 3, the asymptotic forms (^61) and ( ^5|) 



(the different-times generalization of (^)) are exact and the same for all n. In fact, 
asx^l-|-gorg^oo and 7 — equation ( pOD becomes the linear equation (^), 
from which the same powers, but not the amplitudes, can be obtained. 

VI. OTHER SCALING FUNCTIONS IN THE LG APPROACH 

The LG method, implemented in section 5 to evaluate the pair correlation function, 
can be extended to other scaling functions. In this section we evaluate the pressure, 
the average energy density and the energy density correlation function. 

As long as we replace by its saturation form m/|m| (or by m/\m\, for a 
scalar field), the scaling functions will have built in the late-time defect structure. 
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Treating m as a gaussian field, the dynamical dependence of the scaling functions 
is again embodied by 7(r, r/i, 772). In the same spirit as in section 5, we replace 7 
by its large-n limit. In short, we keep in the n-dependence of the scaling properties 
through the gaussian averages over the m vectors, and treat the gaussian moments of 
m in the large-n limit. As mentioned in section 4, the mapping (p^ only determines 
m up to a factor (which may be time-dependent), and thus there is some freedom 
to fix the form of the second moment Sq = (m^). Although the choice So = const, 
would greatly simplify the algebra (e.g. reducing the number of pair contractions 
of gaussian averages containing 0), we find it physically more convenient to regard 
|m| as a length (close to a defect), and thus to keep the scahng form (0), i.e. 
5*0 = 27f/f3. When written in terms of 7, though, the results are independent of 
the choice made. 

The Higgs field energy density (see (|^)) and isotropically averaged pressure are 
given by H 

P{V) = ^2 + M') + (47) 

and scale as like the background p and p. To evaluate p and p within the LG 
approach, we first consider a vector field. In this case, the potential term is negligible 
(and identically zero in the NLSM) and can be ignored. Writing the derivatives of 
(f) in terms of the derivatives of m, expanding the gaussian average into a sum of 
pair contractions, expressing the averages containing Vm0 as the limit 2 ^ 1 of 
derivatives of C(l,2) with respect to the gaussian moments, and treating m in the 
limit n — i> 00, yields 

{^^) = 1) - (ft)') = C,{1, 1) (^^^ (49) 

;(V0t) = ^.(l^l)^ = ^.(l>l)((V0T)^ ' (50) 
where C^(l, 1) = (9C(1, 2)/97(l, 2))2^i = {n - l)/{n - 2) for n > 3, C^(l, 1) = 
ln{L/w) to leading order for n = 2 (where w is the string size, introduced as a 
short-distance cut-off), and 

- 7oo(i,2)2-.l=7oo(i,i)= \ 2 > 



{a - 2)7] 



(V0)2)^ = (ViV2 7oo(l,2))2^i = -V%oo(l,l) = («-l)(<?) , (51) 
where we used (|77|) and Coo = loo- Hence, from (^7D -(pO|), the LG approach gives 
Piv) = ^^^-^^ Pooiv) , P{v) = ^^^-^^ Pooiv) , n>2, (52) 
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p{r]) = \n{L/w) p^{r]) , p{r]) = \ii{L/w) Podv) , n = 2, (53) 

with 

f \ ^ (4 -a) To 1 

In the radiation dominated era, where a = 2, 7oo(l,l) and — V^7oo(l,l) (and poo 
and Poo) have a leading order logarithmic divergence. Their difference, though, is 
finite and gives (T) = Tg/r/^ (see (fTZD ) (and also Poo/Poo = 1/3)- The relevant case, 
however, is the matter dominated era (a = 4), when matter perturbations started 
to grow, yielding p{ri) = with n > 2 and 

Piv) = (6.75 + , n>2. (55) 

V n — 2/ a'^rj'^ 

Although a consistent implementation of the LG method requires the use of the 
NLSM, and thus a vector field, the approach can be extrapolated for a scalar field in 
an elegant manner. This was already done in section 5, where we simply extended the 
results for the scalar field correlation function taking n = 1 in the BPT function (see 
d^)), rather than deriving an equation for C(l,2). The difference for a scalar field 
is that the wall width, w, plays a role in the dynamics, making the scaling functions 
(which contain time-dependent prefactors) differ from their dimensional analysis 
form. Moreover, the potential term in (^71) -(^Sj) has now a relevant contribution 
{V) = (^(V(pj ^ /2a^. A convenient definition for the non-comoving wall width 
(which is constant in time for sharp domain walls) is w = 4/ct, where a is the 
non-comoving (or physical) surface tension given by 

/oo 
dx {d(f)^/dx)^ , (56) 
-oo 

where here (f)w{x) represents a single planar domain wall, and x is a comoving coor- 
dinate normal to the wall. The value of a depends, through 0w(x), on the explicit 
form of the driving potential. In the spirit of the LG approach we exploit the asymp- 
totic mapping 0(m) = sign(m) to perform the gaussian averages, and treat the m- 
correlators in the large-n limit. To evaluate p and p, we write the derivatives of in 
terms of the derivatives of m and expand the gaussian average into a sum of pair con- 
tractions. Noting that dip /dm = cf)' is sharply peaked at m = and that |Vm|m=o = 
1, we get =J^^dmP{m)^'^ = a(rP{0), where P(m) = e-™'/2So/^/2^ is the 
one-point probability distribution for m. Using 0'^ = aa6{m) (which follows from 
([56|)) and integrating by parts, we get ((0'^)") = acr {5"{m)) = aa P"{0). Therefore, 



2^ = a(T^5o/27r 7(1,1) and ((V^)^) = aa^^ So/27r {-V^-f{l,l)), which are the 
analogues of (^S|)-(^np. Since a embodies the extra physical feature of the scalar 
field, we can treat the remaining factors in the large-n limit. Taking 5*0 = ff /Tq 
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(which follows from ([T^) and (^)), 7(i, i) ^ /<f \ and -V^-fil, 1) ^ ((W)^)^, 

\ / o o 

which are then given by (fTTD, we get, from (|I8|)-(P7|), 

/ \ cr 2a — 1 , , , , a 7 — ia , , ^ 

p(r?) = ^== ar]p^{r]) , p r/) = — == — ar] p^{r]) , n = l. 57 

In the radiation dominated era p and p have again a leading logarithmic divergence. 
In the matter dominated era, we obtain 



a 7vC75 1 , , a 3V6775 1 , , 

Piv) = -7^ — z , Piv) = 7k= A ' n = l . 58 

V27r 4 ar] V27r 4 ar] 

We now look to evaluate the correlations between the energy density terms of 
the Higgs field, i.e. and (V0)^. For simplicity we shall restrict to the case of a 
scalar field. Writing the derivatives of in terms of the derivatives of m, expanding 
the gaussian average into a sum of pair contractions, replacing 0'^ by aa6{m), doing 
some gaussian integrals by parts, using (|Il|)-(|T3p, and treating the m-correlators in 
the large-n limit, we obtain (with {XY)^ = (XY) - {X) {Y)) 

(0(1)W)^ + A7n722 = ^^3^ 
{r (7ii-7i)(722-72) + 7^((7i72)^-r7ii722) + 2 (1712+77172)^} , (59) 

(V0(l))2(V0(2))n +AA17A27 ^ 



r(Ai7-7.')(A27-7.')+y(7.'-rAi7A27)+2(r7„+77.')' + ^ , (60) 



(1-72)5/2 

2 , 2r7r 



where A = a"^ a{r]i)r]ia{r]2)r]2/i27rTo), T = (1 - 7^), and 7 = 7oo(l,2), 71 = 
7oo(i,2), 712 = 7oo(i,2), = d-foo/dr, ^„ = d^-foo/dr'^, and 711 = 7oo(i,i), 
Ai7 = — V^7oo(l,l), are given by (|75| ) and (ff^). We have checked that, as ex- 
pected, the results are independent of whether we take 5*0 to be constant or given 
by (|l3). The scaling functions corresponding to ( ^91) and (|60D , normalized in the 
form {XY)^/ {X) (Y), have been plotted in figures 5 and 6, respectively, for the 
matter era. Details and comments are given in the next section and in the figure 
captions. 

VII. SUMMARY AND DISCUSSION 

Two distinct gaussian approaches for the 0{n) Higgs field dynamics, in a fiat ex- 
panding universe, were proposed to evaluate the pair correlation function, and other 
scaling functions. Both theories are based on a non-linear mapping between the 
order parameter 0(r, r/) and an auxiliary field m(r, r]), which varies smoothly in the 



17 



vicinity of the field defects. For simplicity and mathematical convenience, m(r,?7) 
is assumed to be a gaussian random field, yielding an approximate closed scheme to 
evaluate the scaling functions. The field itself, which is effectively discontinuous 
near the defects, is not suitable to be treated as gaussian. 

In the 'soft' field theory of section 3, based on the equation of motion (|^), we 
have followed Mazenko's gaussian approach |]T^ [1^ for model A dynamics, where 



the mapping is defined by the equilibrium profile equation = dVjdcl). In this 

case, m(r, if) is identified as a position vector relative to the nearest field defect. The 
mapping (]T3|), however, is incompatible with the late-time field oscillations in the 
bulk (|l^) (which are absent in purely relaxational systems). By studying the linear 
dynamics of the gaussian moment Co(l, 2), given by equation (|l5l), we can prove that 
this theory is inconsistent, and therefore we have not looked to solve numerically 



the rather complicated equations (|2^) or (^31) . The fact that, despite this intrinsic 
inconsistency, the pair correlation function displays correct physical features, such 
as (p^ and (|26|) , is not, however, a merit of the approximation used. The small-x 
Porod's regime for a scalar field follows from the use of the BPT function (pO|), which 
has built in the late-time defect structure, and the asymptotic power-decay, which 
occurs in the linear (small-/(x)) regime, is universal for all 0{n) 'soft' and 'hard' 
field models. 

In section 4, we have developed a more consistent theory, based on the NLSM, or 
'hard' field, dynamics (p8[). We do not expect the field bulk oscillations ( ]T0| ) to have 
a relevant effect on the scaling properties, so we consistently fix the field magnitude 
to eliminate the previous mapping incompatibility. Also, since the field now evolves 
on the vacuum manifold, the dynamics are independent of the driving potential. 
The auxiliary field is now defined by 0(m) = m/\m\. Although it can still have the 
same interpretation as in section 3, m is only determined up to a factor and we are 
free to choose <m^>. The relation C(7,?t,), between the pair correlation function 
and the normalized m-correlator, is given by the BPT function (^) for a gaussian 
m, and 7 obeys an approximate equation, derivable from (|30D -(p3D. 

Rather then solving this complicated equation for 7, in section 5 we propose 
a fully analytical scheme to evaluate C(l,2). Recognizing that the BPT function 
captures the essential late-time defect structure, we approximate the asymptotic 
field dynamics even further replacing 7 by 700, the exact solution for the limit 
n 00. The pair correlation function is then given (in a symbolic notation) by 
C{1,2)lg = SPT(7oo(l,2),n). Although the NLSM only holds for vector fields, 
the LG approach can be extended to = 1, since it only depends on the large- 
n dynamics, which is the same for both equations (^) and (pS]). In this case the 
scaling properties are evaluated using the mapping (f){m) = sign(m) and the gaussian 
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assumption. The pair correlation function, for instance, is again given by the BPT 
function, with n = 1, and by the same argument we replace 7 by 700. The scaling 
form /lg(3^) Q')) for n = 1, ..,4,oo and in the radiation and matter dominated eras, is 
plotted in: figures 1 and 2, respectively, with different fixed values of q and abscissa 
Xg = 2r/(?7i + 772) = 2x/(l + q), and in figures 3 and 4, respectively, with different 
fixed values of Xs and abscissa q. The normalization is as follows: /(0, 1) = 1 for 
g = 1; for g 7^ 1 we used the time-dependent condition /(x, g)/7oo(a;, g) — 1 as 
X 0, such that curves with different n cut the origin at the same point. 

The LG approach for the Higgs model is the analogue of the Ohta-Jasnow- 
Kawasaki approximate scheme in model A dynamics [0]. In that case, foo{x, 1) = 
exp(— The greater complexity of (^)-(|4l|) is due to the causal condition @ 
which these obey. The main physical features are preserved in this approach: the 
threshold power-law behaviour (|45|) (imposed by causality) is exact, and for n = 1 
a linear Porod's regime, (|43|) , is obtained for a = 4. For a = 2, however, we obtain 
a logarithmic modified Porod's regime, f{x,l) = 1 + 0{x Inx) (slightly apparent 
in figure 1), which is probably an artifact of the LG approach and has no physical 
meaning. This logarithmic correction is absent in the small-x expansion (|24D of 
section 3. 

We have seen in section 5 that all the exact and the gaussian expressions have 
the same limit as n 00. In fact, the gaussian approach becomes exact (for random 
gaussian initial conditions) in this limit since the equation for m becomes linear. This 



equation is derived from the linearised equation for (replacing (f) by m/^/nSo), and 
its form depends on the choice made for Sq. Also, the two gaussian approaches, for 
the 'soft' field and for the NLSM, become equivalent (and exact) as n — > 00 and 
the LG approach could be implemented equally well using either. We find, however, 
that the NLSM provides a more systematic and self-consistent framework for this 
purpose, while the 'soft' field model yields the physical motivation to employ the 
NLSM (and proves useful in the LG calculation of other scaling functions with 
n = 1). 

In section 6 we have extended the LG approach to evaluate other scaling prop- 
erties of the Higgs field. For these cases we do not know how to built closed ap- 
proximate equations like those of sections 3 and 4, and the method proves especially 
useful. If we restricted ourselves to the gaussian approach we could express other 
scaling functions in terms of 7 and its derivatives, but we could not solve for these 
derivatives numerically. For a scalar field we can still use the asymptotic mapping 
(h(m) = m/lml, but we have to account for the non-trivial role of the wall width, 



w, which is inversely proportional to a, the surface tension ([56|). The LG results 
(|5^), (p3D, ([57|), give the average field energy density (^) and pressure ( ^8]) as being 



19 



proportional to Poo? the energy density in the limit n ^ oo. The factor of propor- 
tionality is n-dependent, and is also time-dependent for n = 2 and n = 1. Since 
Poo and Poo, given by (0), have a leading logarithmic divergence at a = 2, we have 
discarded the radiation dominated era, which is a less relevant case in the formation 
of cosmic structure, and next summarize the LG results in the matter era (a = 4). 
With n > 2, (|55| ) gives p = 6.75(1 + l/{n — 2))/a?7f, which compared with the fit 
to simulation results [^: 6.75(1 + (20/9)/(n — 2))/a'^rf shows a fair agreement up 
to a factor ~ 2 in the the correction term. With n = 1, ( pH]) gives p = a const /a?], 
yielding energy density fluctuations growing linearly with time t, rather than having 
a constant value as in the vector case. This well known result means that walls, if 
present, would rapidly dominate the energy of the universe. With > 2 we obtain 
zero pressure, as expected. With n = 1 we get a negative pressure p = — 3p/7, 
yielding a source term p + 3p < 0, which can be regarded as indicating an effective 
domain wall repulsion ||^ in the scaling regime, and is a reflection of domain growth. 
We recall that for an isolated equilibrium domain-wall perpendicular to the x di- 
rection (for which (f)"^ = = d(j)/dy = d(j)/dz), the fleld pressure components along 
each axis are p^ = and Py = Pz = • The pole of (0) at n = 2 is built in the 
approach through the use of a unit vector (i.e. the defect core-size w 0) since, in 
fact, the 'sharpness' of the string cores leads to a logarithmic cut-off given in ( p3D 

m. 



Finally, we have done a LG calculation of the correlations between the energy 
terms (f) and (V0)^, which are the sources for the perturbations in the cosmological 
matter distribution. For simplicity we have restricted ourselves to the case of a scalar 
fleld. In contrast to the vector case, 0^ cannot be regarded as the 'centripetal' energy 
due to the fleld wandering around in the 'vacuum manifold'. At late-times, though, 
0^, or (V0)^, vanish everywhere in the bulk regions and thus probe the presence of 
domain walls (where energy is concentrated) . Using (^8]) and (^^ , we have computed 
the scaling function <^0(1)^0(2)^^ / ^0(1)^^^0(2)^^, given by (|59|), in the matter 
era. Figure 5 shows the results with different flxed values of g = 772/771 and abscissa 
Xs = 2x/{l + q))). Remarkably, as x increases from zero there is a dramatic change 
from large positive values to negative values. We interpret these set of plots as giving 
evidence of domain walls dynamics (in a statistical sense): the correlation peak (for 
flxed q) is displaced along the 'distance' axis as the time separation between the two 
points increases (i.e. as q departs from 1). Its amplitude decrease is dictated by 
statistical incoherence as the points move apart, and its displacement, Xg^peak, must 
be proportional to the typical distance traveled by a wall during the time \rii — 1]2\. 
The equal-time {q = 1) divergence of the peak at the origin (i.e. of ^0^^) is an 
artifact of the absence of a short-distance cut-off in •jooix, 1) as r drops below the 
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wall width w. The scaling function {{V(j){l)y{V(j){2)f)J{{V(f){l)y){{V(j){2)Y), 
given by (|60D , is plotted in figure 6 for the matter era. In this case the peak remains 
at the origin, while its amplitude decreases, as q departs from 1. Since both energy 
density terms probe the presence of domain walls, it is not very clear to us why this 
correlation function is so different from the previous one shown in figure 5. It seems 
that its form for x < 0.5 is entirely dictated by its singular prefactor 1/(1 — 7^)^^^, 
which is plotted in figure 7. 

We conclude by discussing some directions for future work. By linearizing the 
full equation of motion (p!TD, with (|T8|) or (p3D, for the correlation function around 
the scaling solution, it should be possible to show that the scaling solution is a 
stable attractor of the dynamics. In particular, the prescaling regime (e.g. correc- 
tions to scaling) can be described using the gaussian closure schemes of sections 
3 or 4. The early-time behaviour, however, is not accessible within the auxiliary 
field methods utilized here, which assume a well-defined defect structure. The de- 
pendence on the initial state of the system before the phase transition is also of 
interest. In the present work, short-range spatial correlations in the initial state are 
considered, appropriate for a disordered system in equilibrium at high temperature. 
In the context of model A dynamics, it has been shown that the asymptotic scaling 
behavior is modified if sufficiently long-ranged power-law correlations are present 
in the initial state [|2^. Such correlations can also be incorporated in the LG ap- 
proach, through a modification of the large-n solution presented in the appendix. 
Finally, the extension of the functions ( |59D and ( |60D for vector fields, which involves 
extensive calculations, can be used to evaluate the correlations between the matter 



distribution perturbations induced by the Higgs field |TT]] and may, therefore, have 
a direct cosmological interest. 
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APPENDIX: the large-n solution of the NLSM 

For simplicity we shall still take 0^ = 1, which differs from the usual normalization, 
0^ = n, used to solve large-n models. 

To leading order as n — oo the NL factor (|29| ) is replaced by its average (over 
initial conditions), T(r, r/) — ^ (T(r, r/)), which one expects to have the scaling form 
(P7|), i.e. (T(r, r/)) = Tq/t]"^, where Tq is a constant to be determined self-consistently. 
The equation of motion (^) becomes linear and has been solved in momentum 
space, with the following initial conditions, at some early time ^^o > after the SSB 



21 



transition (f){r,rio) is a (gaussian) random unit vector in each initial correlation 
volume, i.e. its Fourier components are white noise correlated, <^0k(^o) '^-kl^o)) = 
A; d(f)k{rio) / dr] — > as k ^ 0, to ensure that dcpir^rjo) /di] respects the assumption 
of homogeneity of the early universe on scales above the horizon. The solution 
obtained is 

u = l + a/2 (62) 
To = 3(2a + l)/4, (63) 

where A^, = 2^T{v + 1) and Ju{z) is a Bessel function of the first kind. The second 
linearly independent solution is ruled out since Yu{kvi) — > oo as k — > 0. From (^TJ) 
one finds the structure factor 

Mvi)-<PMV2)) = — 3— — — — — . (64) 



To obtain the pair correlation function we Fourier transform ([6^). Using the nor- 
malization Coo(0,r7,r7) = 1, we have 



Coc{r,r]i,r]2) = (0(0, ?7i) -(^(r, 772)) 



"^'^ (4(r/i)-0-k(r]2)) e^^ V (4(r]i)-0-k(r/i)) 



{.Vimf^ f Ju{kr]i)J^{kTq2\^ir.\, ^^^^^ 



(66) 



A^i J (27r)3 {kr]iY{kT]2Y 

'Ju{y)V 7r5(a,3/2) 



2-2 r (s±i) 

Clearly, from (|65|), Coo(r, 771, r72) has the scaling form foo{x,q), with x = r/rji and 
Q = V2/vi- To evaluate (|65D we write it in the convenient form 



Using the integral representation of Bessel functions 24 , we have 



7r(i+a) V 2 / J^oo yiyqY 

dy ds dt (1 - s^)''"i/2(l - ^2)^^-1/2 cos(ya;) cos(?/s) cos(ygt) 



Ti(l + a) J-00 Jo Jo 
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\ — — f^ds rdt{l-s^)'^{q^~t'^)'^ {6{x+s-t) + 6{x-s+t) +6{-^+s 

(l + a)g"+^Jo Jo 

/_^^^ (1 - _ (a; _ sf)'^9{s + q- x)9{q + x - s) 



e{l + q-x) 



J ds {1 - s'^)"^\q'^ - [x - s] 



^) 2 



(l + a)g°+2 

alternatively, performing the 'self-similar' transform s ^ x — s in (|68|), 
J(a;,g) = -^^^^^ £"|'ci. (1 - (x - - 3^)^% - .)% + .) 

+ ^^"^rf. (1 - (X - - .^)^ , (70) 



foo{x,q) 



(5, A) = {x + q,x — q) , x<l — q 

= {l,x — q) , |1 — g|<x<l + g 

= (1,-1) , x<q-l, (71) 

Differentiating /(x, q) with respect to x we get, using some of the possible integral 
representations for dl/dx, 

(fl'^ dl{x,q) 



N X dx 

ds (x - s)(g2 -{x- s)2)^(l - s^)^ (72) 



^(1 + g-x) f^,, .,n , „^2^^n 



q^^-i- JA 

a-l 



where = (a + 1)5(0;, 3/2). Expression ([7^) follows from differentiating (|69|). 
The form ([73|) , which follows from (^) and the transformation s — * x — s, or from 
integrating (|72D by parts, is convenient for further differentiation with respect to 
X. Finally, (|7^ is the mean of the previous two, and proves useful at equal-times 
{q = 1) where the factor 1/x gets canceled and higher derivatives with respect to x 
become easier to evaluate. 

By construction J(x, q) must be invariant under interchange of times, i.e. /(x, q) = 
I{x/q,l/q). Since it is not explicitly symmetric, a number of integration variable 
changes and other transformations may be performed in I{x,q) and dl{x,q)/dx 
leading to different equivalent integral representations for foe- However, the expres- 
sions given, with three different integration limits (|7Tj) depending on x and q, admit 
no further simplification. Writing the integrand in, say, (|55|), as ({x,q; s)^~ , it is 
easy to see that ({x, q; s) = (1 — s)(l + s)(g + x — s)(g — x + s) is non-negative and 



23 



bounded only in the regions where both |s| < 1 and x — q<s<x + q, which are 
precisely those yield by (^). Hence, since is non-integer, the integral (|69D runs 
over the whole (bounded) region where the integrand is real. As illustrated by the 
small-x expansion (^21), I{x,q) is singular for a = 2. In fact, at each integration 
limit one (or two, if x = 0) of the radicals in ({x,q;s) vanishes and high enough 
derivatives of the integrand or integration limits will diverge. Up to fourth order, 
however, we get finite derivatives of I{x, q), but since each of the radicals in ( |69D can 
only be differentiated twice, one has to transform the integral, e.g. using r — s — > s 
(or integrating by parts) before doing the third and fourth derivatives. Using these 
methods we find, form (|38D , with 6 = a + 1/2, 

Coo(i,2) = - fiTT^ rdsF,,i-5Co,(l,2)) 

Coo(l,2) = - i^^s rdsqF,,2-SCooil,2)] 
rj2 \Nxq^ Ja j 

Coo(i,2) =^[—1- rrfsgF,,i,2-5'C7^(l,2)-5(r7iCoo(i,2)+r72Coo(l,2) 
r]xr]2 \N xq" JA ^ 

aCoo(l,2) 1/1 /-s C^{\,2)\ 

as 



dr rji \N X q^ J A ^ x 

1 1 fB 



where it is implicit that x < 1 + q, and 

= (g2_c,(a;-s)2)(l-s2)("+l)/2(g2_(a,_,)2){a-3)/2 

F,.3 = -(a + l)s(g2_a(x-s)2)(l-s2)("-i)/2(g2_(a;_s)2)("-3)/2 
= (a + l)(x-s)(l-s2)("-i)/2(g2_(a;_s)2)(«-i)/2 

= (a-l)g(x-s)(l-s2)(°+l)/2(g2_(a._3)2)(a-3)/2 
F^,l,2 = (a'- I)g(x-S)(l-S2)(— l)/2(g2_ (a;_s)2)(a-3)/2 _ (7g) 

We also obtain in the limit 2^1, i.e. r ^ and 1]2 r]i, 



Coo(i,i) = o , vCoo(i,i) = o 

Coo(i,i) = \^ = ~v'cui,i) + ^ 

r]i a — 2 rii 



-V'CUIA) = (77) 
Tji a — 2 
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Figure 1: Field pair-correlation scaling function C(l,2) = f{x,q) in the LG approach for 
the radiation dominated era (see (20) and (38)). Abscissa: Xg = 2 r/(?7i + 772)- Each set 
of lines (crossing at a; = 0) is a collapse of the plots for ra = 1, 2, 3, 4, 00 field components 
with a fixed ratio q = i]2/i]i = 1, 1.5,2.0,2.5,3.0 (bottom). Normalization: f{x, 1) = 1 for 
g = 1, and f{x, q) /joo{x, g) — ?■ 1 as a; — ?■ for all g / 1 and ra, i.e. we have replaced (20) by 
C(7oo(a;, q)) = Too (a;, q) F{a, a; c;j^{x, qy)/F{a, a; c;j^{0, qf ) {a = 1/2, c = (ra + l)/2). 
This time-dependent condition assures that the point where each curve cuts the origin 
is the same for all ra. In all plots f{x,q) = for r > i]i -\- r]2 (causality). The modified 
Porod's regime for ra = 1: f{x, 1) = 1 -\- 0{x ln{x)) as a; — ?■ 0, is an artifact of using the 
large-ra solution. 




Figure 2: The same as in figure 1 but for the matter dominated era. We find the usual 
Porod's regime for ra = 1: /(s, 1) = 1 -|- 0{x) as a; — ?■ (see (43)). 
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Figure 3: Field pair-correlation scaling function C(l,2) = f{x,q) in the LG approach 
for the radiation dominated era (see (20) and (38)). Abscissa: q = i]2/i]i. Each set of 
lines (merging as g — ?■ oo) is a collapse of the plots for n = l,2,3,4,oo with a fixed 
Xg ~ 0.0, 0.3, 0.6, 0.9, 1.2 (bottom). The top curve gives the time-decay at a; = in figure 
1. All curves fall off like as g — ?■ oo (see (46)), and are (by symmetry) invariant 

under the change q — ?■ 1/q. The apparent singularity and correlation increase with q, 
between q = 1 and 2, is an artifact of the time-dependent normalization used in figures 
1 and 2. Using a time-independent normalization, which is then ra-dependent, we find 
that all curves are monotonically decreasing as q departs from 1, but curves with different 
values of Xg and n are difficult to distinguish. 




Figure 4: The same as in figure 3 but for the matter dominated era. 
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Figure 5: Source-source (energy density) pair-correlation scaling function {^(^^(l)(^^(2)y / 

{^(^^(1)^{^(^^(2)^ in the LG approach, given by (59), for a scalar field in the matter domi- 
nated era. Abscissa: Xg = 2 r / {r]i -\- ri^). Each plot is for a different fixed ratio q = ri^/rii. 
from q = 1.0 (top) to g = 2.0 (bottom), with steps A = 0.05. There is interesting evidence 
of domain walls dynamics: the correlation peak (for fixed q) moves along the 'distance' 
axis as time-separation between the two points increases (i.e. as q departs from 1). The 
displacement should be proportional to the typical distance traveled by a wall during time 
|??2 — Vi]- The peak amplitude decreases due to statistical incoherence as the points move 
apart. The equal-time (g = 1) divergence at the origin is an artifact from the assumption 
of 'infinitely sharp' walls (w — > 0). 
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Figure 6: The source-source pair-correlation scaling function in the LG approach 
((V(/>(1))2(V(?!>(2))2)^/ ((V(?!>(1))2)((V(?!>(2))2), given by (60), with the same specifications 
as in figure 5. In this case the correlation peak, while decreasing in amplitude, remains at 
the origin. Below x = 0.5 its form seems to be dictated by its singular prefactor, shown 
in figure 7. 




Figure 7: The prefactor 1/(1 - 7^)®/^ in expressions (59) and (60), as a function of x and 
for fixed values of q in the same range as in figures 5 and 6. 
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